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LONG  TERM  FREQUENCY  STABILITY  ANALYSIS  OF  THE 
GPS  NAVSTAR  6  CESIUM  CLOCK 


INTRODUCTION 

As  part  of  the  Navy  support  to  the  NAVSTAR  Global  Positioning  System  (GPS)  clock  develop¬ 
ment  program,  the  Naval  Research  Laboratory  (NRL)  has  continued  research  and  development  of  pre¬ 
cise  time  and  frequency  standards  [11.  This  report  describes  the  on-orbit  performance  evaluation  of  the 
NAVSTAR  6  cesium  clock. 

The  cesium  clock  in  NAVSTAR  6  is  the  fifth  one  to  be  orbited  in  the  GPS  clock  development 
program.  This  cesium  clock  was  built  under  contract  to  the  Navy  by  Frequency  and  Time  Systems 
(FTS).  The  NAVSTAR  6  cesium  clock  was  activated  on  26  April  1980  (Day  number  117,  1980)  and 
has  been  in  continuous  operation  for  more  than  1  year. 

The  FTS  clock  is  a  preproduction  model  (PPM)  which  is  designated  as  PPM-U.  The  preproduc¬ 
tion  scries  of  cesium  beam  frequency  standards  is  an  evolutionary  development  from  the  prototype 
model  orbited  in  the  Navigation  Technology  Satellite  II  (NTS-2).  The  preproduction  cesium  frequency 
standards  built  by  FTS  are  scheduled  to  be  placed  in  the  NAVSTAR  5,  6,  7,  and  8  spacecraft. 

The  NAVSTAR  6  spacecraft  is  also  known  as  SV  9.  The  SV  identification  is  given  as  part  of  the 
navigation  message. 

GPS  DESCRIPTION 

The  NAVSTAR  GPS  is  a  Department  of  Defense  (DoD)  space-based  system  employing  a  constel¬ 
lation  of  satellites  which  broadcast  signals  that  are  synchronized  in  both  time  and  frequency.  Informa¬ 
tion  is  given  in  the  ephemeris  message  which  can  be  combined  with  measurements  from  four  GPS 
satellites  and  used  to  calculate  the  user’s  instantaneous  position  and  velocity  in  all  three  coordinates,  as 
well  as  precise  time  and  frequency.  Signals  from  the  GPS  satellites  can  be  rapidly  acquired  and  pro¬ 
cessed  independently  of  all  other  systems.  The  precise  time  and  frequency  information  can  be  used  to 
provide  a  common  time  grid  for  worldwide  referencing  of  scientific  laboratory  measurements. 

The  NAVSTAR  GPS  system  comprises  four  segments:  ' 

1.  Control  Segment 

2.  Space  Segment 

3.  User  Segment 

4.  Engineering  Segment 

The  Space,  Control,  and  User  segments  of  GPS  will  be  discussed  with  emphasis  on  factors  that 
are  important  to  this  report. 

The  (iPS  Control  Segment  consists  of  a  Master  Control  Station  (MCS),  located  at  Vandenberg, 
CA,  and  four  remote  Monitor  Sites  (MS),  located  at  Vandenberg,  Hawaii,  Alaska,  and  Guam.  These 
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four  stations  track  the  GPS  spacecraft  vehicles  (SV).  Data  from  these  sites  are  collected  at  the  MCS 
and  processed  to  determine  SV  health,  orbits,  and  clock  offsets. 

The  GPS  Space  Segment  currently  consists  of  six  NAVSTAR  SVs,  with  18  satellites  scheduled  to 
be  operational  in  the  1986-87  time  frame.  These  satellites  are  placed  in  (nominal)  12-h  near-circular 
orbits,  with  occasional  orbit  adjust  maneuvers  which  maintain  a  repeating  ground  track  for  each  SV. 
The  configurations  of  the  NAVSTAR  have  been  under  active  study  since  the  original  recommended 
constellation  [2].  Recent  studies  (3]  have  produced  configurations  that  result  in  a  small  improvement 
in  the  GPS  coverage. 

Each  GPS  spacecraft  broadcasts  spread  spectrum  modulated  signals  that  are  precisely  related  to  the 
on-board  clock.  The  spacecraft  navigation  message  [4]  is  also  modulated  onto  the  signal  at  precise 
epochs  which  aid  in  defining  GPS  time. 

The  "precise,"  or  P-code  modulation  is  generated  at  two  frequencies  in  the  L-band;  these  are 
designated  as  the  L\  and  L2  signals.  The  P-code  signals  are  modulated  at  a  rate  of  10.23  Mbps  (million 
bits  per  second).  The  P-code  modulation  provides  the  capability  for  high  precision  time  difference 
measurements,  and  is  resistant  to  electronic  countermeasures  and  multipath  interferences.  The  P-code 
employs  a  very  long  code  that  is  reset  once  per  week.  The  second  code  which  is  designated  as  the 
coarse/acquisition  code,  or  (C/A)  code,  is  modulated  at  1.023  Mbps  and  repeats  every  millisecond.  It 
provides  a  coarse  signal  that  is  a  factor  of  ten  less  precise  than  the  P-code.  The  C/A  signal  is  quadri- 
phase  modulated  with  the  P-code  and  may  be  rapidly  acquired  by  all  users.  Each  GPS  SV  has  an  atomic 
frequency  standard  that  controls  the  broadcast  frequency  of  each  satellite  to  the  same  nominal  value. 
The  use  of  the  spread  spectrum  modulation  and  separate  codes  for  each  GPS  SV  permits  multiple 
access  to  any  of  the  satellites  that  are  above  the  user’s  horizon. 

A  GPS  user  would  be  required  to  have  an  appropriate  antenna,  receiver,  processor,  and  output 
device  to  receive  the  precise  time  and  time  interval  signals  and  perform  a  navigation  solution.  A  fully 
operational  GPS  user  would  select  four  NAVSTARs  from  the  six  to  nine  satellites  that  would  be  avail¬ 
able  in  such  a  fashion  as  to  minimize  the  Geometrical  Dilution  Of  Precision  (GDOP),  a  quantity  [5,21 
that  relates  to  the  navigation  accuracy  available  from  GPS.  The  user  would  acquire  and  lock  the 
receiver  to  signals  broadcast  from  four  of  the  SVs,  and  then  make  simultaneous  measurements  of  time 
difference  (pseudo-range)  and  frequency  difference  (pseudo-range-rate)  between  each  of  the  SV  clocks 
and  the  receiver  clock.  The  user  would  then  use  the  four  pseudo-range  measurements  to  calculate 
clock  offset,  latitude,  longitude,  and  height.  The  four  pseudo-range-rate  measurements  would  be  used 
to  calculate  frequency  offset  and  velocity  in  all  three  components. 

LONG  TERM  FREQUENCY  STABILITY 

The  GPS  system  is  capable  of  providing  instantaneous  precise  navigation  because  the  satellite 
clocks  are  synchronized  in  time  and  frequency.  Therefore  a  fundamental  measure  of  system  perfor¬ 
mance  is  given  by  the  long  term  frequency  stability  of  each  of  the  SV  clocks. 

The  technique  for  analyzing  frequency  stability  performance  of  orbiting  clocks  and  frequency  stan¬ 
dards  was  developed  at  NRL  in  1975  [6]. 

This  procedure  has  evolved  into  an  analytical  procedure  depicted  in  Fig.  1.  Each  msqor  com¬ 
ponent  of  the  long  term  frequency  stability  analyses  is  described  in  this  report. 
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Fig.  1  —  GPS  frequency  stability  analysis  procedure  for  analyzing  on-orbit  clock  performance 


FREQUENCY  STABILITY  MODEL 


The  Allan  variance  was  adopted  by  the  IEEE  as  the  recommended  measure  of  frequency  stability. 
Reference  7  presents  a  theoretical  development  that  results  in  a  relationship  between  the  expected 
value  of  the  standard  deviation  of  the  frequency  fluctuation  for  any  finite  number  of  data  samples  and 
the  infinite  time  average  of  the  standard  deviation.  Equation  (1)  presents  the  Allan  variance  expression 
for  M  frequency  samples  with  sample  period  T  equal  to  the  sampling  time,  r: 


a-y( 2,  r)  — 


J  M-\ 

Tm^Y) 


(vk+i  -  yJ1 

2 


(1) 


The  average  frequency  values  yk,  calculated  from  pairs  of  clock  offsets,  At,  separated  by  sample 
time,  r,  are  given  by: 


yk 


T 


(2) 


The  clock  offset,  satellite  orbit,  and  other  variables  must  be  measured  or  estimated  from  the 
pseudo-range  measurements. 
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CLOCK  DIFFERENCE  MEASUREMENTS 

Pseudo-range  (PR)  and  accumulated  delta  pseudo-range  (ADR)  measurements  are  taken  between 
the  NAVSTAR  SV  clock  and  the  MS  clock  using  a  spread  spectrum  receiver.  The  measurements  are 
taken  once  every  6  s  and  then  aggregated  and  smoothed  once  per  1 '  min.  Figure  2  presents  a  plot  of  a 
typical  pseudo-range  signature  obtained  from  a  single  NAVSTAR  pass  over  a  monitor  station.  Each 
measurement  is  corrected  for  equipment  delay,  ionospheric  delay,  tropospheric  delay,  earth  rotation, 
and  relativistic  effects.  Then  the  data  are  edited  and  smoothed  using  the  predicted  SV  ephemeris  to  cal¬ 
culate  the  geometric  delay.  The  clock  offset  at  the  midpoint  of  the  15-min  data  span  is  estimated  by 
using  both  the  pseudo-range  and  the  pseudo-range  rate  measurements,  which  are  fitted  to  a  cubic  poly¬ 
nomial  with  epoch  at  time  corresponding  to  the  midpoint  of  the  data. 


Fig.  2  —  Typical  NAVSTAR  pseudo-range 
signature,  plotted  every  IS  min 


The  pseudo-range  measurements  are  resolved  to  1/64  of  a  P-code  chip,  which  corresponds  to 
1.5  ns  or  46  cm  in  range.  Nominal  values  for  the  pseudo-range  noise  levels  are  o-PR  =  1.3  m  for  the  L\ 
measurements  and  aPR  -  2.0  m  for  the  L2  measurements.  The  L\  and  L2  measurements  are  com¬ 
bined  to  correct  for  ionospheric  refraction,  which  results  in  an  increase  to  4.53  m  for  the  corrected 
pseudo-range  measurement.  The  accumulated  della  pseudo-range  measurement  noise  levels  are 
0.31  cm  for  L\  and  0.56  cm  for  Lv  These  measurements  are  also  combined  to  correct  for  the  accumu¬ 
lated  pseudo-range  measurements.  The  smoothing  procedure  uses  the  ADR  measurements  to  aid  the 
PR  smoothing  of  each  15  min  segment  of  data.  This  process  results  in  a  smoothed  pseudo-range  meas¬ 
urement  noise  level  of  18.5  cm. 

The  equation  that  relates  the  pseudo-range  measurements  to  the  clock  difference  between  the 
NAVSTAR  SV  and  the  MS  is: 

PR  "  R  +  c(f\j  s  —  fsv)  4-  cA  tA  -he,  (3) 

where 


PR  -  the  measured  pseudo-range, 

R  —  the  slant,  or  geometric  range,  from  the  SV  (at  the  time  of  transmission)  to  the  MS  (at 
the  time  of  reception), 
c  -  the  speed  of  light, 
rMS  -  the  MS  clock  time, 
rSv  —  the  SV  clock  time, 

tA  -  ionospheric,  tropospheric,  and  relativistic  delay,  with  corrections  for  antenna  and  equip¬ 
ment  delays, 
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and 


e  -  the  measurement  error. 

The  clock  difference,  denoted  by  Atk  for  the  kth  measurement  is  obtained  by  rearranging  Eq.  (3)  into 

Atk  =  (rSv  ~  *ms)  =  R/c  +  AtA  +  e/c  -  PR/c.  (4) 

The  particular  evaluation  of  Eq.  (4)  that  will  be  used  in  this  report  is  obtained  by  designating  the 
NAVSTAR  6  by  SV  9  and  by  designating  the  Vandenberg  Monitor  Site  as  VMS.  This  particular  case  of 
Eq.  (4)  is  given  by: 

Atk  =  (SV9  —  VMS).  (5) 


SMOOTHED  ORBIT  ESTIMATION 


All  of  the  smoothed  pseudo-range  measurements  are  collected  from  the  four  GPS  monitor  sites 
for  1  week.  The  Naval  Surface  Weapons  Center  (NSWC)  then  estimates  a  smoothed  NAVSTAR  orbit 
using  an  orbit  estimation  program  that  extensively  models  the  dynamics  of  the  satellite  motion,  includ¬ 
ing  solar  radiation  pressure,  and  orbit  adjust  maneuvers.  The  NSWC  post-fit  ephemeris  calculations 
employ  the  highly  redundant  set  of  range-difference  [5]  values,  which  are  calculated  from  the  smoothed 
15  min  pseudo-range  measurements. 

The  purpose  of  the  smoothed  orbit  estimation  is  to  separate  the  clock  and  orbital  components  by 
modeling  the  clock  as  a  constant,  (but  unknown)  frequency  during  the  1-week  span.  The  model 
includes  an  (unknown)  aging  rate,  which  may  be  used  for  frequency  standards  that  exhibit  aging.  The 
model  also  is  capable  of  segmenting  the  clock  bias  solution  to  allow  for  frequency  adjustments  of  the 
MS  or  SV  clock. 


The  clock  differences  used  for  analyzing  the  spacecraft  clock  incorporate  the  smoothed  orbit  and 
the  set  of  15-min  pseudo-range  measurements  to  calculate  the  clock  difference  at  the  time  of  each 
measurement  according  to  Eq.  (4).  The  clock  differences  for  each  NAVSTAR  pass  are  then  used  to 
estimate  the  clock  differences  at  the  time-of-closest-approach  (TCA)  of  the  NAVSTAR  SV  over  the 
monitor  site.  This  procedure  results  in  either  one  or  two  points  per  day.  The  NAVSTAR  orbit  and  the 
monitor  site  location  determine  whether  one  or  two  points  per  day  will  be  available.  For  NAVSTAR  6, 
one  pass  per  day  is  available  from  the  Vandenberg  Monitor  Station. 


The  evaluation  of  the  clock  difference  at  the  TCA  point  minimizes  the  effect  of  the  NAVSTAR 
orbit  estimation  for  along-the-satellite-track  and  out-of-plane  errors.  However  this  procedure  does  not 
reduce  the  effect  of  radial  orbit  errors.  Hence,  the  estimate  of  radial  orbit  error  will  be  one  of  the  fac¬ 
tors  that  limits  the  accuracy  of  the  long  term  frequency  stability  analysis.  The  effect  of  the  radial  error 
on  the  frequency  stability  is  given  by: 


Vyir) 


'Siv  RK 
cr 


(6) 


where 


<trr  -  standard  deviation  estimate  of  radial  component  of  orbit  error, 
c  -  the  speed  of  light, 

r  —  sample  time, 

and  the  units  of  <rRR,  c,  and  t  result  in  a  dimensionless  number  for  <ry(-r). 
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NAVSTAR  6  RESULTS 

The  clock  differences  between  the  NAVSTAR  6  cesium  clock  and  the  VMS  clock  are  presented  in 
Fig.  3.  Each  "X"  symbol  corresponds  to  a  single  measurement  obtained  from  the  smoothed  15-min 
pseudo-range  measurement.  Twenty-three  points  are  plotted  for  this  NAVSTAR  6  pass  over  the  VMS. 
These  23  points  are  analyzed,  and  a  subset  of  these  data  are  used  to  estimate  the  clock  offset  at  the 
TCA  point. 


Fig.  3  —  Single  Pass  Clock  difference  between  NAVSTAR-6  (SV9)  and 
the  GPS  Vandenberg  Monitor  Site  (VMS) 


j  The  clock  difference,  which  is  denoted  by  (SV9  -  VMS),  corresponds  to  starting  a  time  interval 

measurement  with  the  NAVSTAR  6  clock,  and  stopping  with  the  VMS  clock  (with  corrections  for  the 
orbit  and  other  delays).  The  clock  offset  changed  by  approximately  30  ns  during  the  SV  pass.  The 
clock  offset  presented  here  represents  a  clock  difference  that  may  be  processed  further  to  produce 
"GPS"  time.  The  slope  of  these  measurements  indicates  a  fractional  frequency  offset  of —1.08  x  10“ 13 
fractional  between  the  NAVSTAR  6  clock  and  the  VMS  clock.  The  magnitude  of  this  offset  is  normally 
|  what  would  be  expected  [8]  after  the  correction  for  the  relativistic  clock  effect. 

The  clock  differences  for  one  week  are  presented  by  Fig.  4.  In  this  figure,  each  "X"  symbol 
denotes  one  clock  difference  obtained  from  the  smoothed  15-min  pseudo-range  measurements.  There 
are  seven  groups  of  "X"  symbols,  each  one  corresponding  to  a  single  NAVSTAR  6  pass  that  was 
observed  by  the  VMS.  The  slope  of  the  clock  differences  for  this  I-week  segment  is  -1.56  x  10~12. 

The  clock  differences  for  the  entire  100-day  data  span  are  presented  in  Fig.  5.  Each  vertical  mark 
corresponds  to  the  clock  difference  evaluated  at  the  TCA  point  of  a  NAVSTAR  pass  over  the  VMS. 

Figure  5  indicates  a  total  change  in  clock  difference  of  about  10  ps  in  100  days,  or  approximately 
0.1  fts/day,  or  equivalently,  a  fractional  frequency  offset  of  -1.16  x  10-12.  Close  analysis  of  Fig.  5  also 
indicates  that  a  small  frequency  change,  on  the  order  of  5  x  10-13,  occurred  between  days  120  and  150. 
This  small  frequency  shift  will  be  further  analyzed  before  computing  the  frequency  stability. 

The  frequency  differences  for  sample  times  of  1,  3,  and  10  days  are  presented  in  Figs.  6,  7,  and  8. 
Analysis  of  the  results  indicates  that  the  total  frequency  change  between  the  two  cesium  standards  dur¬ 
ing  the  100-day  span  was  on  the  order  of  1  x  10-IJ.  Figure  7  indicates  three  small  frequency  changes 
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Fig.  4  —  One  week  of  clock  differences  between 
NAVSTAR-6  (SV9)  and  the  GPS  Vandenberg 
Monitor  Site  (VMS) 
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Fig.  5  —  100  days  of  clock  differences  between 
NAVSTAR-6  (SV9)  and  the  GPS  Vandenberg  Moni¬ 
tor  Site  (VMS) 
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Fig.  7  —  100  days  of  frequency  differences  between 
NAVSTAR-6  (SV9)  and  GPS  Vandenberg  Monitor 
Site  (VMS)  for  a  3-day  sample  time 


Fig.  8—100  days  of  frequency  differences  between 
NAVSTAR-6  (SV9)  and  the  GPS  Vandenberg  Moni¬ 
tor  Site  (VMS)  for  a  10-day  sample  time 


on  the  order  of  3  x  10-13  occurred,  with  the  majority  of  the  frequency  differences  on  the  order  of  a  few 
parts  in  10(14).  Analysis  of  other  GPS  data  (not  included  in  this  report)  indicates  that  the  VMS 
cesium  clock  was  responsible  for  the  largest  frequency  changes  analyzed. 

The  frequency  measurements  were  then  used  to  calculate  the  Allan  variance  for  sample  times 
varying  from  1  to  10  days.  An  interleaving  (also  called  overlapping)  data  processing  technique  was 
used  in  order  to  obtain  maximal  use  of  the  data.  For  instance,  with  a  sample  time  of  1  day  and  the  set 
of  clock  differences  [M\,  At2,  At,,  A/4}  two  variances  were  calculated.  The  first  variance  used  the  sub¬ 
set  {A/|,  At 2,  Ar3}  and  the  second  variance  was  calculated  using  the  subset  {Af2,  At),  Af4).  Thus,  the 
two  a-ylr )  values  have  the  subset  {Ar2.  A/3}  in  common. 

Frequency  stability  estimates,  for  the  combined  clock  and  cphemeris,  are  presented  in  Fig.  9.  The 
1 -sigma  upper  and  lower  confidence  limits  for  v/t )  have  been  calculated  assuming  flicker  noise  FM 
for  the  noise  process.  The  results  are  summarized  by  Table  1. 

These  results  may  be  characterized  by  four  segments,  using  a  model  given  by: 

(TyW)  «■  a  rM, 


(7) 


where  the  coefficient  a  and  the  exponent  p  are  experimentally  determined  from  the  data  by  use  of  a 
model  described  in  [71. 

The  coefficients  and  sample  times  for  the  combined  clock  and  ephemeris  frequency  stability  are 
summarized  in  Table  2. 


Table  2  —  Coefficients  of  the  Time  Domain 
Model  for  the  Allan  Variance 

Sample  Time  Coefficient  "a"  Exponent  V" 


(days) 

1  <  r  <  3 
3  <  t  <  5 
5  <  t  <  9 
9  <  r  <  10 


PP  10(13) 
1.10 
0.90 
0.33 
1.30 


-0.18 

0.00 

0.63 

0.00 
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Comparison  of  the  exponents  with  the  type  of  noise  process  identifiable  in  [7]  in  atomic  frequency 
standards  indicates  that  the  segment  for  sample  times  from  3  to  5  days  is  classified  as  flicker  noise  fre¬ 
quency  modulation  (FM).  The  segment  for  sample  times  of  9  to  10  days  may  also  be  classified  as 
flicker  noise  FM.  The  other  two  segments  can  not  readily  be  classified,  however  the  first  segment  has 
an  exponent  of  -0.18  which  is  close  to  flicker  noise  FM,  and  the  third  segment  has  a  slope  of  0.633, 
which  is  close  to  that  expected  from  a  random  walk  in  frequency  (with  an  exponent  of  1.0). 

The  error  sources  that  are  believed  to  be  most  significant  in  limiting  this  analysis  are: 

1.  The  use  of  a  single  monitor  station  frequency  standard, 

2.  the  radial  component  of  orbit  error. 

The  first  factor  can  be  reduced  by  incorporating  multiple  frequency  standards  into  the  analysis  at  a 
single  monitor  station,  or  by  analyzing  the  spacecraft  clock  performance  using  multiple  monitor  sites. 
For  example,  the  use  of  the  other  three  GPS  monitor  sites  will  permit  the  identification  of  frequency 
changes  such  as  evidenced  in  Fig.  7.  A  time  scale  could  then  be  formed  in  a  manner  similar  to  that 
described  in  [8],  The  effect  of  the  radial  orbit  error  can  be  estimated  by  using  Eq.  (6).  For  the  best  fit 
ephemerides  used  in  this  analysis,  the  average  radial  error  was  2.1  m,  which  corresponds  to  1.4  x  10~13 
for  a  sample  time  of  1  day.  Additional  orbit  smoothing  could  produce  better  estimates  for  the  orbit; 
however,  the  - 1  value  of  the  r  exponent  in  Eq.  (6)  is  unchanged  by  additional  smoothing.  Ultimately, 
the  18-cm  noise  level  obtained  from  the  smoothed  15-min  measurements  yields  a  lower  limit  of 
1.2  x  10“ 14  for  a  1-day  sample  time. 

CONCLUSIONS 

1.  The  frequency  stability  performance  of  the  NAVSTAR  6  cesium  clock,  to  date,  is  accept¬ 
able  and  within  specifications. 

2.  The  combined  (spacecraft  clock,  single  monitor  station  clock,  and  ephemeris)  frequency 
stability  is  equal  to,  or  less  than,  1.3  x  10~u  for  sample  times  of  1  to  10  days. 
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